Human conjunctiva organoids to study ocular surface homeostasis and disease.

The conjunctival epithelium covering the eye contains two main cell types: mucus-producing goblet cells and water-secreting keratinocytes, which present mucins on their apical surface. Here, we describe long-term expanding organoids and air-liquid interface representing mouse and human conjunctiva. A single-cell RNA expression atlas of primary and cultured human conjunctiva reveals that keratinocytes express multiple antimicrobial peptides and identifies conjunctival tuft cells. IL-4/-13 exposure increases goblet and tuft cell differentiation and drastically modifies the conjunctiva secretome. Human NGFR+ basal cells are identified as bipotent conjunctiva stem cells. Conjunctival cultures can be infected by herpes simplex virus 1 (HSV1), human adenovirus 8 (hAdV8), and SARS-CoV-2. HSV1 infection was reversed by acyclovir addition, whereas hAdV8 infection, which lacks an approved drug therapy, was inhibited by cidofovir. We document transcriptional programs induced by HSV1 and hAdV8. Finally, conjunctival organoids can be transplanted. Together, human conjunctiva organoid cultures enable the study of conjunctival (patho)-physiology.

Regulation of lymphoid-myeloid lineage bias through regnase-1/3-mediated control of Nfkbiz.

ABSTRACT: Regulation of lineage biases in hematopoietic stem and progenitor cells (HSPCs) is pivotal for balanced hematopoietic output. However, little is known about the mechanism behind lineage choice in HSPCs. Here, we show that messenger RNA (mRNA) decay factors regnase-1 (Reg1; Zc3h12a) and regnase-3 (Reg3; Zc3h12c) are essential for determining lymphoid fate and restricting myeloid differentiation in HSPCs. Loss of Reg1 and Reg3 resulted in severe impairment of lymphopoiesis and a mild increase in myelopoiesis in the bone marrow. Single-cell RNA sequencing analysis revealed that Reg1 and Reg3 regulate lineage directions in HSPCs via the control of a set of myeloid-related genes. Reg1- and Reg3-mediated control of mRNA encoding Nfkbiz, a transcriptional and epigenetic regulator, was essential for balancing lymphoid/myeloid lineage output in HSPCs in vivo. Furthermore, single-cell assay for transposase-accessible chromatin sequencing analysis revealed that Reg1 and Reg3 control the epigenetic landscape on myeloid-related gene loci in early stage HSPCs via Nfkbiz. Consistently, an antisense oligonucleotide designed to inhibit Reg1- and Reg3-mediated Nfkbiz mRNA degradation primed hematopoietic stem cells toward myeloid lineages by enhancing Nfkbiz expression. Collectively, the collaboration between posttranscriptional control and chromatin remodeling by the Reg1/Reg3-Nfkbiz axis governs HSPC lineage biases, ultimately dictating the fate of lymphoid vs myeloid differentiation.

Spatial and Temporal Mapping of Breast Cancer Lung Metastases Identify TREM2 Macrophages as Regulators of the Metastatic Boundary.

Cancer mortality primarily stems from metastatic recurrence, emphasizing the urgent need for developing effective metastasis-targeted immunotherapies. To better understand the cellular and molecular events shaping metastatic niches, we used a spontaneous breast cancer lung metastasis model to create a single-cell atlas spanning different metastatic stages and regions. We found that premetastatic lungs are infiltrated by inflammatory neutrophils and monocytes, followed by the accumulation of suppressive macrophages with the emergence of metastases. Spatial profiling revealed that metastasis-associated immune cells were present in the metastasis core, with the exception of TREM2+ regulatory macrophages uniquely enriched at the metastatic invasive margin, consistent across both murine models and human patient samples. These regulatory macrophages (Mreg) contribute to the formation of an immune-suppressive niche, cloaking tumor cells from immune surveillance. Our study provides a compendium of immune cell dynamics across metastatic stages and niches, informing the development of metastasis-targeting immunotherapies. SIGNIFICANCE: Temporal and spatial single-cell analysis of metastasis stages revealed new players in modulating immune surveillance and suppression. Our study highlights distinct populations of TREM2 macrophages as modulators of the microenvironment in metastasis, and as the key immune determinant defining metastatic niches, pointing to myeloid checkpoints to improve therapeutic strategies. This article is featured in Selected Articles from This Issue, p. 2489.

The transcriptional and regulatory identity of erythropoietin producing cells.

Erythropoietin (Epo) is the master regulator of erythropoiesis and oxygen homeostasis. Despite its physiological importance, the molecular and genomic contexts of the cells responsible for renal Epo production remain unclear, limiting more-effective therapies for anemia. Here, we performed single-cell RNA and transposase-accessible chromatin (ATAC) sequencing of an Epo reporter mouse to molecularly identify Epo-producing cells under hypoxic conditions. Our data indicate that a distinct population of kidney stroma, which we term Norn cells, is the major source of endocrine Epo production in mice. We use these datasets to identify the markers, signaling pathways and transcriptional circuits characteristic of Norn cells. Using single-cell RNA sequencing and RNA in situ hybridization in human kidney tissues, we further provide evidence that this cell population is conserved in humans. These preliminary findings open new avenues to functionally dissect EPO gene regulation in health and disease and may serve as groundwork to improve erythropoiesis-stimulating therapies.

A single-cell census of mouse limb development identifies complex spatiotemporal dynamics of skeleton formation.

Limb development has long served as a model system for coordinated spatial patterning of progenitor cells. Here, we identify a population of naive limb progenitors and show that they differentiate progressively to form the skeleton in a complex, non-consecutive, three-dimensional pattern. Single-cell RNA sequencing of the developing mouse forelimb identified three progenitor states: naive, proximal, and autopodial, as well as Msx1 as a marker for the naive progenitors. In vivo lineage tracing confirmed this role and localized the naive progenitors to the outer margin of the limb, along the anterior-posterior axis. Sequential pulse-chase experiments showed that the progressive transition of Msx1+ naive progenitors into proximal and autopodial progenitors coincides with their differentiation to Sox9+ chondroprogenitors, which occurs along all the forming skeletal segments. Indeed, tracking the spatiotemporal sequence of differentiation showed that the skeleton forms progressively in a complex pattern. These findings suggest an alternative model for limb skeleton development.

Cellular and metabolic characteristics of pre-leukemic hematopoietic progenitors with GATA2 haploinsufficiency.

Mono-allelic germline disruptions of the transcription factor GATA2 result in a propensity for developing myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML), affecting more than 85% of carriers. How a partial loss of GATA2 functionality enables leukemic transformation years later is unclear. This question has remained unsolved mainly due to the lack of informative models, as Gata2 heterozygote mice do not develop hematologic malignancies. Here we show that two different germline Gata2 mutations (TgErg/Gata2het and TgErg/Gata2L359V) accelerate AML in mice expressing the human hematopoietic stem cell regulator ERG. Analysis of Erg/Gata2het fetal liver and bone marrow-derived hematopoietic cells revealed a distinct pre-leukemic phenotype. This was characterized by enhanced transition from stem to progenitor state, increased proliferation, and a striking mitochondrial phenotype, consisting of highly expressed oxidative-phosphorylation-related gene sets, elevated oxygen consumption rates, and notably, markedly distorted mitochondrial morphology. Importantly, the same mitochondrial gene-expression signature was observed in human AML harboring GATA2 aberrations. Similar to the observations in mice, non-leukemic bone marrows from children with germline GATA2 mutation demonstrated marked mitochondrial abnormalities. Thus, we observed the tumor suppressive effects of GATA2 in two germline Gata2 genetic mouse models. As oncogenic mutations often accumulate with age, GATA2 deficiency-mediated priming of hematopoietic cells for oncogenic transformation may explain the earlier occurrence of MDS/AML in patients with GATA2 germline mutation. The mitochondrial phenotype is a potential therapeutic opportunity for the prevention of leukemic transformation in these patients.

LGR5 expressing skin fibroblasts define a major cellular hub perturbed in scleroderma.

Systemic sclerosis (scleroderma, SSc) is an incurable autoimmune disease with high morbidity and mortality rates. Here, we conducted a population-scale single-cell genomic analysis of skin and blood samples of 56 healthy controls and 97 SSc patients at different stages of the disease. We found immune compartment dysfunction only in a specific subtype of diffuse SSc patients but global dysregulation of the stromal compartment, particularly in a previously undefined subset of LGR5+-scleroderma-associated fibroblasts (ScAFs). ScAFs are perturbed morphologically and molecularly in SSc patients. Single-cell multiome profiling of stromal cells revealed ScAF-specific markers, pathways, regulatory elements, and transcription factors underlining disease development. Systematic analysis of these molecular features with clinical metadata associates specific ScAF targets with disease pathogenesis and SSc clinical traits. Our high-resolution atlas of the sclerodermatous skin spectrum will enable a paradigm shift in the understanding of SSc disease and facilitate the development of biomarkers and therapeutic strategies.

The interaction of CD4+ helper T cells with dendritic cells shapes the tumor microenvironment and immune checkpoint blockade response.

Despite their key regulatory role and therapeutic potency, the molecular signatures of interactions between T cells and antigen-presenting myeloid cells within the tumor microenvironment remain poorly characterized. Here, we systematically characterize these interactions using RNA sequencing of physically interacting cells (PIC-seq) and find that CD4+PD-1+CXCL13+ T cells are a major interacting hub with antigen-presenting cells in the tumor microenvironment of human non-small cell lung carcinoma. We define this clonally expanded, tumor-specific and conserved T-cell subset as T-helper tumor (Tht) cells. Reconstitution of Tht cells in vitro and in an ovalbumin-specific αß TCR CD4+ T-cell mouse model, shows that the Tht program is primed in tumor-draining lymph nodes by dendritic cells presenting tumor antigens, and that their function is important for harnessing the antitumor response of anti-PD-1 treatment. Our molecular and functional findings support the modulation of Tht-dendritic cell interaction checkpoints as a major interventional strategy in immunotherapy.

Physically interacting beta-delta pairs in the regenerating pancreas revealed by single-cell sequencing.

OBJECTIVES: Until recently, communication between neighboring cells in islets of Langerhans was overlooked by genomic technologies, which require rigorous tissue dissociation into single cells. METHODS: We utilize sorting of physically interacting cells (PICs) with single-cell RNA-sequencing to systematically map cellular interactions in the endocrine pancreas after pancreatectomy. RESULTS: The pancreas cellular landscape features pancreatectomy associated heterogeneity of beta-cells, including an interaction-specific program between paired beta and delta-cells. CONCLUSIONS: Our analysis suggests that the particular cluster of beta-cells that pairs with delta-cells benefits from stress protection, implying that the interaction between beta- and delta-cells might safeguard against pancreatectomy associated challenges. The work encourages testing the potential relevance of physically-interacting beta-delta-cells also in diabetes mellitus.

COVID-19 in Patients with Hematologic Malignancies: Clinical Manifestations, Persistence, and Immune Response.

BACKGROUND: The clinical presentation of coronavirus disease 19 (COVID-19) is the result of intricate interactions between the novel coronavirus and the immune system. In patients with hematologic malignancies (HM), these interactions dramatically change the clinical course and outcomes of COVID-19. SUMMARY: Patients with HM and COVID-19 are at an increased risk for prolonged viral shedding, more protracted and severe presentation, and death, even when compared to other immunocompromised hosts. HM (e.g., multiple myeloma, chronic lymphocytic leukemia) and anticancer treatments (e.g., anti-CD20 agents) that impair humoral immunity markedly increase the risk of severe COVID-19 as well as protracted viral shedding and possibly longer infectivity. Cytokine release syndrome (CRS) is an important player in the pathophysiology of severe and fatal COVID-19. Treatments targeting specific cytokines involved in CRS such as interleukin-6 and Janus kinase have proven beneficial in COVID-19 patients but were not assessed specifically in HM patients. Although neutropenia (as well as neutrophilia) was associated with increased COVID-19 mortality, granulocyte colony-stimulating factors were not beneficial in patients with COVID-19 and may have been associated with worse outcomes. Decreased levels of T lymphocytes and especially decreased CD4+ counts, and depletion of CD8+ lymphocytes, are a hallmark of severe COVID-19, and even more so among patients with HM, underlying the important role of T-helper dysfunction in severe COVID-19. In HM patients with intact cellular immunity, robust T-cell responses may compensate for an impaired humoral immune system. Further prospective studies are needed to evaluate the mechanisms of severe COVID-19 among patients with HM and assess the efficacy of new immunomodulating COVID-19 treatments in this population. KEY MESSAGES: Understanding the immunopathology of COVID-19 has greatly benefited from the previous research in patients with HM. So far, no COVID-19 treatments were properly evaluated in patients with HM. Patients with HM should be included in future RCTs assessing treatments for COVID-19.

XCR1+ type 1 conventional dendritic cells drive liver pathology in non-alcoholic steatohepatitis.

Non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH) are prevalent liver conditions that underlie the development of life-threatening cirrhosis, liver failure and liver cancer. Chronic necro-inflammation is a critical factor in development of NASH, yet the cellular and molecular mechanisms of immune dysregulation in this disease are poorly understood. Here, using single-cell transcriptomic analysis, we comprehensively profiled the immune composition of the mouse liver during NASH. We identified a significant pathology-associated increase in hepatic conventional dendritic cells (cDCs) and further defined their source as NASH-induced boost in cycling of cDC progenitors in the bone marrow. Analysis of blood and liver from patients on the NAFLD/NASH spectrum showed that type 1 cDCs (cDC1) were more abundant and activated in disease. Sequencing of physically interacting cDC-T cell pairs from liver-draining lymph nodes revealed that cDCs in NASH promote inflammatory T cell reprogramming, previously associated with NASH worsening. Finally, depletion of cDC1 in XCR1DTA mice or using anti-XCL1-blocking antibody attenuated liver pathology in NASH mouse models. Overall, our study provides a comprehensive characterization of cDC biology in NASH and identifies XCR1+ cDC1 as an important driver of liver pathology.

Meningeal lymphoid structures are activated under acute and chronic spinal cord pathologies.

Tertiary lymphoid structures (TLS) are organized aggregates of B and T cells formed ectopically during different stages of life in response to inflammation, infection, or cancer. Here, we describe formation of structures reminiscent of TLS in the spinal cord meninges under several central nervous system (CNS) pathologies. After acute spinal cord injury, B and T lymphocytes locally aggregate within the meninges to form TLS-like structures, and continue to accumulate during the late phase of the response to the injury, with a negative impact on subsequent pathological conditions, such as experimental autoimmune encephalomyelitis. Using a chronic model of spinal cord pathology, the mSOD1 mouse model of amyotrophic lateral sclerosis, we further showed by single-cell RNA-sequencing that a meningeal lymphocyte niche forms, with a unique organization and activation state, including accumulation of pre-B cells in the spinal cord meninges. Such a response was not found in the CNS-draining cervical lymph nodes. The present findings suggest that a special immune response develops in the meninges during various neurological pathologies in the CNS, a possible reflection of its immune privileged nature.

Publisher Correction: Cxcl10+ monocytes define a pathogenic subset in the central nervous system during autoimmune neuroinflammation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Author Correction: Cxcl10+ monocytes define a pathogenic subset in the central nervous system during autoimmune neuroinflammation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Host-Viral Infection Maps Reveal Signatures of Severe COVID-19 Patients.

Viruses are a constant threat to global health as highlighted by the current COVID-19 pandemic. Currently, lack of data underlying how the human host interacts with viruses, including the SARS-CoV-2 virus, limits effective therapeutic intervention. We introduce Viral-Track, a computational method that globally scans unmapped single-cell RNA sequencing (scRNA-seq) data for the presence of viral RNA, enabling transcriptional cell sorting of infected versus bystander cells. We demonstrate the sensitivity and specificity of Viral-Track to systematically detect viruses from multiple models of infection, including hepatitis B virus, in an unsupervised manner. Applying Viral-Track to bronchoalveloar-lavage samples from severe and mild COVID-19 patients reveals a dramatic impact of the virus on the immune system of severe patients compared to mild cases. Viral-Track detects an unexpected co-infection of the human metapneumovirus, present mainly in monocytes perturbed in type-I interferon (IFN)-signaling. Viral-Track provides a robust technology for dissecting the mechanisms of viral-infection and pathology.

Cxcl10+ monocytes define a pathogenic subset in the central nervous system during autoimmune neuroinflammation.

Multiple sclerosis (MS) is characterized by pathological inflammation that results from the recruitment of lymphoid and myeloid immune cells from the blood into the brain. Due to subset heterogeneity, defining the functional roles of the various cell subsets in acute and chronic stages of MS has been challenging. Here, we used index and transcriptional single-cell sorting to characterize the mononuclear phagocytes that infiltrate the central nervous system from the periphery in mice with experimentally induced autoimmune encephalomyelitis, a model of MS. We identified eight monocyte and three dendritic cell subsets at acute and chronic disease stages in which the defined transcriptional programs pointed toward distinct functions. Monocyte-specific cell ablation identified Cxcl10+ and Saa3+ monocytic subsets with a pathogenic potential. Transfer experiments with different monocyte and precursor subsets indicated that these Cxcl10+ and Saa3+ pathogenic cells were not derived from Ly6C+ monocytes but from early myeloid cell progenitors. These results suggest that blocking specific pathogenic monocytic subsets, including Cxcl10+ and Saa3+ monocytes, could be used for targeted therapeutic interventions.

Dissecting cellular crosstalk by sequencing physically interacting cells.

Crosstalk between neighboring cells underlies many biological processes, including cell signaling, proliferation and differentiation. Current single-cell genomic technologies profile each cell separately after tissue dissociation, losing information on cell-cell interactions. In the present study, we present an approach for sequencing physically interacting cells (PIC-seq), which combines cell sorting of physically interacting cells (PICs) with single-cell RNA-sequencing. Using computational modeling, PIC-seq systematically maps in situ cellular interactions and characterizes their molecular crosstalk. We apply PIC-seq to interrogate diverse interactions including immune-epithelial PICs in neonatal murine lungs. Focusing on interactions between T cells and dendritic cells (DCs) in vitro and in vivo, we map T cell-DC interaction preferences, and discover regulatory T cells as a major T cell subtype interacting with DCs in mouse draining lymph nodes. Analysis of T cell-DC pairs reveals an interaction-specific program between pathogen-presenting migratory DCs and T cells. PIC-seq provides a direct and broadly applicable technology to characterize intercellular interaction-specific pathways at high resolution.